Low-energy spectra in the t-J type models at light doping 
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Abstract 

Based on a variational approach, we propose that there are two kinds of low 
energy states in the t-J type models at low doping. In a quasi-particle state an 
unpaired spin bound to a hole with a well-defined momentum can be excited 
with spin waves. The resulting state shows a suppression of antiferromagnetic 
order around the hole with the profile of a spin bag. These spin-bag states with 
spin and charge or hole separated form a continuum of low-energy excitations. 
Very different properties predicted by these two kinds of states explain a 
number of anomalous results observed in the exact diagonalization studies on 
small clusters up to 32 sites. 



Typeset using REVTgX 



Soon after the discovery of the high T c superconducting cuprates, the t-J model was 
proposed [1] as the prototype to examine the phenomena. Since then it has been shown in 
theories [2,3] and experiments [4,5] that extra hoppings beyond nearest neighbors (n.n.) are 
also important to the low-energy features of the cuprates. Enormous theoretical effort has 
been devoted to predict the low-energy spectra of these t-J type models. But the strong cor- 
relation associated with these models has made perturbative approaches ineffective. With 
the many different mean-field theories presented so far, there is little consensus to the "cor- 
rect" description of low-energy spectra of these models. To sort out the proper theory it 
is important to first have a detail comparison of the predictions with the results of exact 
calculations. 

Recently, significant progress has been made at very low doping. Several different nu- 
merical techniques, such as exact diagonalization (ED) studies [6-9], self-consistent Born 
approximation (SCBA) [10] and the Green function Monte Carlo method [11], all have ob- 
tained similar results for the energy-momentum dispersion relation of a single hole doped 
into t-J type models. The results agree fairly well with experiments [5]. In addition, a 
mean-field or variational wave function (VWF) [12] constructed from the half-filled Mott in- 
sulating state with antiferromagnetic long-range order (AF LRO) has also obtained a similar 
success. In a recent paper [3] this single-hole VWF is generalized to treat systems with mul- 
tiple holes or electrons. The new set of VWF's easily explained angle-resolved photoemission 
spectroscopy (ARPES) results [4] for Fermi pockets around (tc/2,tc/2) and (vr,0) for hole- 
doped and electron-doped systems, respectively. It [3] also reproduced unusual patterns in 
momentum distribution functions (MDF's) of the ground state calculated by the ED method 
[13,14] for one and two holes in 32 sites. For these VWF's, doped holes or electrons behave 
like quasi-particles (QP's). In this paper we will show that in addition to these QP states, 
there is a continuum of charge excitations described by spin-charge separated states in the 
spectra of t-J type models. The presence of the two kinds of states is carefully examined by 
explaining several anomalous results reported by ED studies for clusters up to 32 sites. 

As shown by Lee et al. [3], the ground state in the presence of a few doped holes 
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or electrons could be described by a VWF constructed from the half-filled Mott insulat- 
ing state. At half-filling, the ground state of the Heisenberg Hamiltonian is described 
fairly accurately by including three mean-field paramters [12,16]: the staggered magneti- 
zation m s =(S%)=—{S^), where the lattice is divided into A and B sublattices; the uni- 
form bond order parameters x=(J2a c L c jo-)j an d rf-wave resonating valence bond (<i-RVB) 
order A=(c^Qf — CjfC^) if % and j are n.n. sites in x-direction and —A for ^-direction. 
Without d-EVB order, the mean-field Hamiltonian has lower and upper spin-density-wave 
(SDW) bands with operators aka=«kCk CT + cr/3 k Ck+Q CT , and &ka=-C/3kCk CT + ^kCk+Qa, re- 
spectively Here Q—(n, 7r), a k =|[l — (e k /£ k )] and /3k=|[l+(ek/£k)]- Energy dispersions 
for the two SDW bands are ±£ k = ±(ek + {Jfn s ) 2 )^ with e k =— jJx(cosk x + cosk y ). In- 
clusion of d-RVB pairing for electrons on total N s sites, the VWF for the ground state 
has the form \^o)=Pe[Ek(A k a^al kl + B k b^bl kl )] N ^ 2 \0), where A k =(E k + &)/A k and 
B\ t =~(E ii — £ k )/A k with _Ek=(^k + A k )^ 2 an d the constraint of one electron per site en- 
foced by P e . Here A k =| JAd k with d^coskj. — cosk r Notice that the sum in \^f ) is taken 
over the sublattice Brillouin zone (SBZ). 

In the presence of doped holes or electrons, we consider the t-J type model including 
longer-ranged hoppings, with amplitues t' for the 2nd n.n. and t" the 3rd n.n.. By applying 
a particle- hole transformation [3,6] we can treat hole- and electron-doped cases in the same 
manner. However, here we will just concentrate on the hole-doped cases with J/t—0.3 = 
-t'/t and t"/t=0.2. 

When a hole is doped or an electron is removed from \^o), a P arr must be broken with 
an unpaired spin left. Thus it is quite natural to have the following VWF [3,12] for a single 
doped hole, e.g., with a lone up spin 

\^ l (q s ))=P d cl 1 

[ E (^a^ + i?^^)]^/^-^), 

[k^q h ] 

where the hole momentum is excluded from the sum if q/j is within the SBZ, otherwise, 
qt — Q is excluded. Pd here enforces the constraint of no doubly occupied sites. When we 
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choose the unpaired-spin momentum q s to be either the same as the hole momentum q/j 
or q^+Q, this VWF is equivalent to the Lee-Shih [12] wave function. Variational energies 
calculated vary with q/j [3,12]. The energy dispersions for t-J and t-t'-t"-J models are plotted 
as filled circles in Fig. 1(a) and (b), respectively. For both models, the ground state with 
one hole has momentum (ir/2,ir/2) = Q/2. As shown in Ref. [3], these dispersion relations 
are still followed when hole number is increased. The holes in these wave functions behave 
just like QP's, hence we denote |\l/i(q/ l =q s )) = |\l/? p ). 

There are only two variational parameters: A/x and m s /x in our VWF's. The extended 
hoppings t' and t" are not used as variational parameters in both \^f) and The effect 

of t is included in the RVB uniform bond x- Since t' and t" are compatible with AF LRO, 
there is little effect for them to be included in our VWF. 

Clearly, the choice of unpaired spin to have the same momentum as the hole, i.e. q s =q/i, 
is a special case for \%>i). If we choose q s ^ q^, then not only the electron pair at q^ and 
— q/j is excluded in the sum in the pair at q s and — q s is also affected. Hence we expect 
it to be higher in energy. To make a distinction from the afore-mentioned QP states l^f), 
these states will be denoted as the spin-bag (SB) states, i.e. \^\(q s 7^ q&)) = \^f)- The 
variational energies as a function of q s for |\Pf b ) with q h =Q/2 are plotted as empty circles 
in Figs. 1(a) and (b) for the t-J and t-t'-t"-J models, respectively. Many SB states could be 
constructed with the same q s but different q^. While it is possible to have the SB states of 
even lower energy with q ft =(37r/4, 0) in the t-J model (gray circles in Fig. 1(a)), they are 
of higher energies than that of SB states with <\ h ={n/2, ±7r/2) in the t-t'-t"-J case. There 
are many such states, actually an infinite number of them for an infinite system, forming a 
continuum, as schematically illustrated by the shaded regions in Fig. 1(a) and (b). 

There is an intuitive way to understand the difference between SB and QP states. The 
spin excitations of the QP states can be easily constructed by applying the spin operators, 



jection operator P^. In the linear spin-wave theory k' is the momentum of the spin wave. 
The particular term inlcuded in the sum of q' with q' equal to the momentum of the un- 
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Notice that these operators commute with the pro- 



paired spin q s =qh changes the QP state \^f) to the SB ones \^f(q s =q h + k')). Thus 
the SB states are actually just spin- wave excitations of the QP state with the same hole 
mementum. In Fig. 1(c) the difference of variational energies between the SB states and 
QP state with q/j=Q/2 for the t-t'-t"-J model is plotted as a function of the difference of 
momentum k'=q s — q h . The dotted line is the prediction of energy-momentum dispersion 
relation of linear spin-wave theory [17]. The slight differences between the two results at 
k'=Q/2 and k'=(7r, 0) are due to the hole-renormalization effect [18]. SB states represent 
spin excitations of the QP states. The empty/gray circles in Figs. 1(a) and (b) are the 
lowest spin excitation energies of the ground state [19]. 

The two VWF's also give very different spin and hole correlations. When the hole is at 
the A sublattice, the correlation is SH A (r)=J2 ie A(-l) i+r (^ll^iSi^!) / EieA(^i\ n i\^i}, 
where is the hole- number operator at site i; similarly for the case when hole is at B 
sublattice. In our convention, the up spin prefers to be on sublattice A. If the hole and 
unpaired spin are uncorrelated, SHa,b equals to the value of uniform staggered magnetiza- 
tion. Values of SHa for the QP state (filled triangles) with q s =q/ l =(0,0) and the SB state 
(empty triangles) with q s =(0,0) and q^=Q/2 [20], respectively, are plotted as a function 
of distance from the hole in Fig. 2. Both states have an unpaired down spin and a single 
hole. The spin configurations around the hole are clearly different for the two states. The 
spin magnetization right next to the hole in the QP state \^f) has values larger than the 
uniform background, 0.368. Thus the unpaired spin is bound to the hole. On the other 
hand, for the SB state the magnetization is suppressed around the hole, this is similar to 
the idea of a SB first proposed in Ref. [21]. The unpaired spin bound to the hole in the 
QP state is here being excited by the spin-wave excitation and becomes unbound in the SB 
state. This may be viwed as a spin-charge separated state. 

The spin-charge separation observed in the SB states has many interesting consequences. 
In ED studies, it has been found [15] that the lowest energy state at (n, 0) for the t- 
t'-t"-J model is very different from that of the t-J model. The spin-spin correlation 
across the hole changes from ferromagnetic (FM) to AF when t' is turned on. This re- 
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suit can now be understood as that the lowest energy state at (it, 0) for the t-J model is 
the QP state \^f(q s =q h =(n, 0))) as shown in Fig. 1(a), but it changes to a SB state 
(q s =(7r,0),q ft =Q/2)) for the t-t'-t"-J model as shown in Fig. 1(b). In Table I we list 
correlations obtained between pairs of spins around the doped hole for QP and SB states. 
The correlation is defined as C s ,s>(q s ) = E i (^T ,sb ((l s )\n'lS i+ s ■ S i+l5 ,|*f ,sb (q s )) with 5 and 5' 
denoting two sites around the hole [15]. While the lower-energy QP states at q s =q/ l =(7r, 0) 
and (0,0) behave as expected for a system with AF LRO, i.e. FM for spins at the same 
sublattice (pairs a and b) and AF otherwise, the result for a SB state with q s =(vr, 0) and 
q?i=Q/2 shows AF correlation at the same sublattice This is exactly the behavior observed 
in ED results [15]. Our result in Fig. 1 shows that at (0,0) the lowest energy state remains 
to be the QP state even when t' and t" are included. The spin correlations are thus not 
changed, this is also consistent with what is found by Tohyama et al. [15]. All these results 
clearly support our identification of low-energy spectra to be consisted of QP and SB states. 

TABLE I. C$,6' (q s ) calculated for an 8x8 lattice using different VWF's with momenta 
indicated. The first and third rows are for QP states while others are for SB ones, a to 
e are different pairs of sites defined in the inset of Fig. 2. Positive (negative) values mean 
FM(AF) correlations. 



q* 


q^ 


a b c d e 


(0,0) 


(0,0) 


0.188 0.188 0.202 -0.273 -0.264 


(0,0) 


Q/2 


-0.0288 -0.0254 -0.0302 -0.203 -0.195 


M) 


(tt,0) 


0.123 0.15 0.071 -0.353 -0.279 


M) 


Q/2 


-0.0313 -0.0085 -0.002 -0.1921 -0.212 



Another important difference between QP and SB states is in their MDF (n a (k)). The 
results can actually be predicted. Since the hole momentum q h is excluded from the VWF 
we naturally expect (n-j-(k)) to have a smaller value or a dip at k=q/j and k=q^ + Q 
than its neighbors, similarly for (ri|(k)) at k=— q^ and k=— q fe + Q. However, in a QP state 
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with an up spin at momentum q s =qh, then (nf(k=q^)) is increased and there is no more 
a dip. As an example, the MDF obtained by the QP state with q/ l =q s =(7r, 0) is listed in 
Fig. 3(a). Because of the symmetry, only results for one quadrant of the BZ are shown. 
At each k the upper(lower) number is for up(down) spin. Now for a SB state with a lone 
up spin at q s (^ q/J, the original spin at q^ in the QP state is excited and placed at q s , 
then (n-|-(k=q/i)) still has a dip. Fig. 3(b) shows the MDF obtained for the SB state with 
q s =(7r,0) and q ft =Q/2 [20]. Results are here for 32 sites and (A/%, m s /x)=(0.1, 0.05). 

This behavior of the MDF's is indeed found in the exact results of the t-t'-t"-J and 
t-J models on 32 sites. MDF's obtained by the ED method for the lowest energy state at 
(7r, 0) of the t-J and t-t'-t"-J models are shown in Fig. 3(c) and (d), respectively. The nice 
qualitative agreement achieved between 3(a) and (c) as well as between 3(b) and (d) [22] 
re-affirms our results shown in Fig. 1(a) and (b): the lowest energy state at (ir, 0) for a 
single hole is a QP state for the t-J model and a SB state for the t-t'-t"-J model. 

Another consequence of this switch from a QP state to a SB one is the drastic change 
of the spectral weight, ^ k =|(^i(k)|c k(J |^o)| 2 /(^o|cLcka|^o)- Using |*f (q=(?r, 0))) we ob- 
tained Z k = 0.475, it vanishes when we use the SB state |^f 6 (q s =(7r, 0), q h =Q/2)). This 
is consistent with exact results for the t-J model (Z k =0.34) [7] and t-t'-t"-J model (Z k =0) 
[23]. In addition, spectral weights of the lowest energy states of both models at (ir,7i) and 
(37r/4, 37r/4) are either exactly zero or very small. This is consistent with our identification 
that states at both momenta are SB ones. Since c ka \^f ), unlike the SB state, has momenta 
of the hole and unpaired spin related, it has a negligible overlap with the SB state. By 
contrast, states at (n/2,n/2) and (tt/A,tc/A) remain to be QP states in both t-J and t-t'- 
t"-J models, hence large spectral weights are expected. It is noted that our QP (SB) states 
predict much larger (small) spectral weights in comparison with that of the exact 32 sites. 
This discrepency is partly due to the fact that we have AF LRO in our VWF's while total 
spin is a good quantum number in the exact results. Another reason is that due to the 
projection operator P d our QP states and SB states with the same quantum numbers (total 
momentum and total S z ) are actually not orthogonal to each other although they have very 
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small overlap. But there are many SB states in the continuum that could couple with a 
particular QP state. Hence, when the QP state has energies very close to the continuum, 
the spectral weight of the QP state is diluted by the coupling with SB states. This effect 
makes the quantitative prediction of spectral weight difficult. We leave this issue for the 
future work. 

In summary, based on a mean-field theory with AF and rf-RVB order parameters we have 
proposed that at low doping there are two kinds of low energy states for t-J type models. 
The single-hole QP states have a well-defined energy dispersion. By exciting the QP states 
with spin waves we obtain a continuum of SB states. The unpaired spin is separated from 
the hole in the SB states. A number of physical properties predicted by these two kinds of 
states are in good agreement with the exact results obtained by ED studies. Although our 
emphasis in this paper is to show the solid theoretical support of these two kinds of states, 
there are also experimental evidences. In Ref. [3], QP states were shown to explain well the 
single-hole dispersion observed by ARPES. However the overall variation pattern of spectral 
weights observed in the ARPES experiment on Ca 2 Cu0 2 Cl2 [4] is naturally understood with 
the presence of SB states: notable lowest energy peaks are only observed in small regions of 
k-space, e.g. near (tc/2,tt/2) and (7r/2, 0) where QP states have lower energy than the SB 
states (Fig. 1(b)). More comparison with experiments is in progress. 

We are grateful to Profs. N. Nagaosa and C.-T. Shih for invaluable discussions and their 
supporting data. TKL is supported by the grant NSC91-2112-M-001-011 (R.O.C.). PWL is 
supported by the Hong Kong RGC grant number HKUST6159/01P. 
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FIGURES 

FIG. 1. Variational energies calculated for the (a) t-J and (b) t-t'-t"-J model Hamilto- 
nians for one hole on an 8x8 lattice by applying our TWF's. Filled circles, connected by 
solid lines, are VMC results using the |^^ p (q s )} discussed in the text; empty(gray) circles by 
|*! 6 (cUi q/i = Q/2((37r/4, 0))). Shaded regions indicate a possible continuum for an infinite system. 
A minus sign has been multiplied to all the data shown here, (c) Difference of variational energies 
between the QP ground state at Q/2 and the SB states at q s =Q/2 + k' in (b) as a function of k'. 
The dotted line is the prediction of energy dispersion of linear spin- wave theory [17]. Results here 
are obtained with parameters (A/%, m s /x)=(0. 25, 0.125). 

FIG. 2. Spin-hole correlation functions, SHa (defined in the text), for an 8x8 lattice calculated 
using the QP state (filled triangles) with q s =q/j=(0,0) and the SB state (empty triangles) with 
q s =(0, 0) and q^=Q/2 [20]. Inset: pairs of sites, denoted by letters a to e, where spin-spin 
correlations listed in Table I are computed. Parameters are the same as used in Fig. 1. 

FIG. 3. Momentum distribution functions in one quadrant of k space the BZ for 32 sites 
obtained by (a) |*f (q=(vr, 0))), (b) (q s = (vr, 0), q ft =Q/2)), (c) ED results for (it, 0) state of 
the t-J model and (d) t-t'-t"-J model. The upper number is result for up spin and lower for down 
spin. 
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